ORIGINAL ARTICLE 



Folia Forestalia PolonUa, series A, 2013, Vol. 55 (2), 58-71 



DOI: 10.2478/ffp-2013-0007 



Improved methods of classification of multispectral aerial 
photographs: evaluation of floodplain forests in the 
inundation area of the Danube 

Tomas Bucha 1 M, Martin Slavik 2 

1 National Forest Centre, Forest Research Institute, T.G. Masaryka 22, 960 92 Zvolen, Slovakia, e-mail: bucha@nlcsk.org 

2 Czech University of Life Sciences, Faculty of Forestry and Wood Sciences, Kamycka 129, CZ - 165 21 Praha 6 - Suchdol, 
Czech Republic 



Abstract 

The Gabcikovo hydroelectric power plant has significantly influenced Danube water regime, thus the condition of 
floodplain forests in the region. Forest condition has been regularly monitored since 1995 using aerial photos. The 
subject of this study was to improve the procedure of floodplain forest health evaluation based on digital multispec- 
tral aerial images. Firstly, the forest mask was created with overall accuracy 89%, and next, tree health was evaluated 
using defoliation as health indicator. We applied orthogonal transformation of 4 original bands of multispectral im- 
agery into two-dimensional space. Marginal values of digital numbers (DN) of the first component (New Synthetic 
Channel - NSCl) were defined by fully foliated willow and poplar. The second component (NSC2) was optimised 
for damage estimation. Calculated DN values of NSC2 represented a perpendicular distance from the line of DN 
values of the first component. The distance from the line was proportionate to tree damage extent in a given pixel. 
We generated linear regression model between pair values ofNSC2 and defoliation evaluated for 38 trees in the field, 
respectively, from aerial photos. A decline prediction resulted in r-square equal 0.86. Finally, we used the model to 
predict defoliation for each picture element (pixel) of the component NSC2. 
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Introduction 

The Waterworks Gabcikovo was put into operation in 
1992. Its construction ensured energy exploitation of 
the Danube; waterways transport at the critical part of 
the ford and safety against floods. A crucial Danube part 
- from Hrusov to Sap - was re-routed into a 31 km deri- 



vational canal. Negative impacts on floodplain forests 
are connected with this solution caused by groundwa- 
ter level lowering due to drainage effect of the Danube 
riverbed. Therefore, suitable environment for existence 
of floodplain forests towards preservation of their pro- 
ductive potential was created by implementing hydro- 
technical modifications. Forest health state is monitored 
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in transboundary cooperation following the agreement 
between the governments of the Slovak Republic and 
Hungary that was concluded in 1995. 

Aerial photos are commonly used for ecological 
monitoring (Morgan et al. 2010). For the purposes of 
evaluation of forest condition based on infrared pho- 
tographs, visual interpretation keys were developed 
(Gross 2000; Ciesla 2000). Damage evaluation is based 
on a combination of shape, texture, and colour char- 
acteristics. Automated procedures of the classification 
were developed for digital aerial multispectral images. 
There have been carried out assessments concerning 
forest damage classification based on the use of the 
maximum likelihood method (Meddens et al. 2011), re- 
gression analysis (Bucha et al. 2009), neural networks 
(Klobucar 2010), non-parametric k-nearest neighbour 
method (Eigirdas et al. 2013). In automated approaches, 
spectral differences between healthy and damaged trees 
are used for identifying damage extent. High absorption 
in the blue and red parts of the spectrum and moderately 
increased reflectivity in the green part of the spectrum 
are typical for healthy vegetation. Chronic damage 
causes deterioration of chloroplasts, and this change 
in tree physiology residts in leaves turning yellow, and 
then the maximum reflectivity moves from the green 
part to the red part of the spectrum (Zarco-Tejada et al. 
2001). High and stable reflectivity in the near infrared 
range of the spectrum is typical for healthy vegeta- 
tion. Our measurements of poplar leaves showed high 
reflectivity for dry and brown foliage as well. Unlike 
healthy vegetation it was not uniform but was gradu- 
ally increasing in a range from 700 to 1100 nm, simi- 
larly to the pattern observed by Ahern (1988). Slaton et 
al. (2001) proved relation between infrared reflectance 
and changes of leaf intra-cellular and cellular struc- 
ture. Studies on vegetation response to stress showed 
that within wavelength range from 400 to 850 nm the 
maximum difference in reflectance between control 
and stressed vegetation was close to 700 nm (Carter and 
Knapp 2001; Masaitis et al. 2013). 

Delineation of tree crowns is necessary for a forest 
damage classification. Many approaches have been so 
far proposed to separate tree crowns. Pitkanen (2001), 
Pouliot at al. (2002), Sumbera (2003), Wang et al. (2006) 
proposed the automated delineation of tree crowns for 
black and white, infrared and multispectral aerial im- 
ages. The authors assumed that the top of the tree can 



be detected at the highest brightness value. The edge 
of the crown, on the contrary, represents the minimum 
brightness value. 

Hirschmugl et al. (2007) used stereo images to de- 
rive Digital Model of Terrain (DMT). This was used 
to distinguish individual tree crowns and combined 
with procedures based on image spectral properties, 
namely the maximum and minimum reflectance on 
the top and edge of the crown. Wolf and Heipke (2007) 
proposed a procedure for automated delineation of tree 
crowns from infrared images combined with DMT- 
based fuzzy logic and approximation of tree crown by 
ellipsis. Bijker et al. (2010) modelled the probability 
profiles of single tree crowns with Gaussian functions, 
resulting in delineation of tree crown objects. Jing et 
al. (2012) proposed a method for individual tree crown 
delineation based on multi-scale filtering and imagery 
segmentation. Pontius et al. (2008) used supervised 
classification - Spectral Angle Mapper method - to 
identify non-forested areas. 

Differences of radiation reflectivity depend on tree 
species. In the Danube inundation area, there prevail 
poplars and willows. Higher reflectivity in the visible 
infrared part of the spectrum is typical for willows and 
white poplars when compared to that of black poplar. 
This makes the model for damage estimation harder 
to derive. The problem was solved by application of 
multiple regression analysis in previous cycles of for- 
est monitoring in the Danube inundation. Apart from 
aerial photographs, there were also inserted into the 
model as a dependent variable, data on a representation 
of tree species taken from the forest management plan 
(FMP) (Rasi and Bucha 2001) or damage classification 
was carried out separately for willow and poplar forests 
(Bucha et al. 2009). 

The aim of this study was to improve the procedure 
of floodplain forests health evaluation based on com- 
mercially available aerial photographs. Specifically, our 
objectives were to: 

1. delineate the part of tree crowns on which damage 
would be evaluated; 

2. objectify digital classification of forest condition by 
transformation of original bands into the compo- 
nent optimised for damage estimation; 

3. use the component (derived in the objective 2) in 
conjunction with ground control data to predict for- 
est health condition over large contiguous areas. 
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Material and methods 
S tudy area 

The evaluated area of the Danube left bank inundation 
is located between the original Danube riverbed and 
the canal of the Gabcikovo Waterworks (fig. 1). The 
area is characteristic of developed river branch system, 
vast complexes of floodplain forests and alluvial mead- 
ows. Forest area is approximately 3000 ha. According 
to FMP (2005), floodplain forest species composition 
within the area observed is as follows: 62% cultivated 
poplars (Populus sp.), 10% domestic poplars ( Populus 



nigra L. a P. alba L.), 17% willows ( Salix alba, S. fra- 
gilis L.), 6% ash trees (especially Fraxinus angustifolia 
VAF1L.), and 5% other deciduous trees, especially Eng- 
lish oaks ( Quercus robur L.), acacias ( Robinia pseudo- 
acacia L.), sycamore maples ( Acer pseudoplatanus L.) 
and alders (Alnus sp.). 

Aerial photography 

The basic parameters of the aerial photography are 
described in tab. 1. Organisational and administrative 
works, as well as all the legal requisites of the aerial 
imaging, were ensured by the company Photomap Ltd. 
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Fig. 1. Study area 
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Kosice. Photography was performed on 10 lh Septem- 
ber 2011 by the company ARGUS GEO SYSTEM, Ltd. 
Hradec Kralove using a digital camera Vexcel Ultra- 
CamX. 



Tab. 1 . Selected parameters for aerial photography in 20 1 1 



Locality 


Hrusov - Sap 


Scale 


1 : 41 800 


Endlap 


60% 


Sidelap 


No - photography in 1 strip 


Date of Photography 


1 0 th September 2011 


Time of Photography 


11:00-15:00 


Type of camera 


UltraCamX - digital multispectral 
camera 


Spatial resolution 


30 x 30 cm 


Radiometric resolution 


Blue - Green - Red - Near- 
infrared band 


Total imaging area 


ca. 20 000 ha 



The photographs were supplied in tif format in 8 bit 
resolution, in the 3 rd level of processing, which means 
that the multispectral channels (blue, green, red and 
near-infrared) were spatially modified from 90 x 90 cm 
into 30 x 30 cm resolution according to the panchro- 
matic band. The area of Hrusov - Sap was depicted in 
30 aerial photographs obtained with overlap 60%. All 
the photographs were ortho-rectified using Image sta- 
tion software. The total root-mean square error (RMSE) 
did not exceed the value of ± 1 m at individual photo- 
graphs after their transformation to the national S-JTSK 
coordinate system. All 30 photographs were used for 
creation of stereo-pairs for visual assessment of tree 
defoliation for the purposes of derivation of regressive 
model for damage estimation. The 13 photographs, spe- 
cifically the images 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 
26 and 28 (i.e. every second photo from Slovak territo- 
ry) were used for the classification of forests condition. 

Tree defoliation estimation 

Defoliation (loss of assimilation organs) was used as the 
main indicator and basic visual symptom of tree health. 
The indicator expresses percentage ratio of missing as- 
similation organs to the ideal status of the assimilation 
apparatus, i.e. to a sample fully foliaged tree in given 
conditions. 



The aim of tree defoliation assessment was to ob- 
tain support information for the evaluation of the aerial 
photographs. Defoliation was assessed during field sur- 
vey and by visual interpretation of the photos on the 
monitor screen. 

During field survey there were chosen 11 areas with 
healthy to severely damaged trees. The assessment of 
defoliation was performed on 1 14 trees by one evaluator 
at the end of August 2011, i.e. before image acquisition. 
Defoliation was assessed using binocular - in percents 
with rounding off to 5%, following the methodology 
of the international programme UN-ECE ICP Forests 
(2006). Only 6 from 114 trees were unambiguously 
identified on the aerial photos. These trees served as 
a calibration sample for visual interpretation of defolia- 
tion from photo stereo-pairs as well as were included 
into the regression model for damage estimation. De- 
foliation of further 32 trees was assessed based on ste- 
reo-pairs and methodology described in EU’s manual 
(Gross 2000) by the evaluator who carried out a terres- 
trial assessment. In total, 38 trees were finally used in 
deriving the regression model for estimation of forest 
damage based on the aerial photos: 15 willows, 15 pop- 
lars and 8 dead trees. 

Spectral measurements 

Spectral reflectivity measurements were conducted 
with Licor-1800 spectroradiometer. External integra- 
tion sphere LI-1800-12 was used for the measure- 
ments in a spectral range 400-1100 nm with a 2 nm 
step. Leaves of black poplar, white poplar, grey poplar, 
white willow, quaking aspen and dry bark of aspen 
were analyzed. The leaves were collected in the ar- 
boretum Borova hora in Zvolen from the lower part 
of tree crown in the first decade of August. Measure- 
ments were carried out during app. 60 minutes after 
sample collection. The aim of the measurements was 
to compare determined spectral reflectivity with DN 
values from the aerial photographs and to use the in- 
formation obtained for verification of the proposed 
method of forest damage classification. 

Creation of the forest mask based on aerial 
photographs 

Distinguishing forest from other land cover categories, 
that are overall called forestless, is a significant opera- 
tion because all the characteristics of forest health are 
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Fig. 2. Working flow of the forest mask creation and tree-crown defoliation estimation 



determined based only on pixels classified as forest. We 
created the forest mask for each photograph. Two basic 
processes of classification: object-based approach and 
pixel-based approach were combined (fig. 2). 

Image segmentation: object-based classification 
approach 

Object oriented approach was based on image segmen- 
tation -iirffiTiiore homogeneous objects and performed 
using ^Cognition software. After that we carried out 
combined automated and manual classification of seg- 
ments into two classes: non-forest (forestless segments) 
and other (segments with forest). 



The aim of image segmentation was to distinguish 
forest from other land cover classes. Prior to segmen- 
tation we reduced the image resolution from 30 cm to 
2.4 m which made the process of picture segmentation 
faster (about 30 seconds for 2115 x 2035 pixels). Seg- 
mentation parameters were set empirically based on 
quality of visual examination of the results from several 
image segmentations. We chose the parameters as fol- 
lows: shape criterion = 0.1, compactness criterion = 0.5 
and size criterion = 100. 

A low value of the shape criterion means that colour 
has a great importance at the objects creation. The com- 
pactness criterion defined integrity of a resultant object. 
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The value 0.5 (middle of the interval 0-1) was selected 
based on factual occurrence of various object shapes 
- from compact to non-compact. The size criterion de- 
termined the size of the resultant object. The selected 
value of 100 resulted to number of objects from 390 to 
670 in individual photographs. This range allowed ef- 
fective manual classification of various forestless cat- 
egories. The “line selection” eCognition functionality 
which enables selection of several objects with an only 
draw of a cursor was used. This method of object se- 
lection was used particularly for the classification of 
heterogeneous segments as village with gardens, agri- 
culture areas with hedgerows and individual trees and 
non-arable areas with shrub succession. Homogenous 
segments with crops with similar spectral properties as 
forests as well as water bodies outside forest complexes 
were also selected by this method and then classified 
manually into the forestless category. 

Automated classification was used for water bodies 
and shadows occurring spatially inside forest complex- 
es. The size of these segments is usually smaller what 
makes the manual classification slow. So water bodies 
as shadows are distinctly spectrally different from for- 
est when their DN value in infrared band is lower then 
90. The segment with DN ired < 90 was assigned to the 
forestless category. 

Distinction of tree-crown: pixel-based classification 
approach 

Subsequently, we distinguished tree crowns from other 
categories (shadows, small gaps in vegetation, roads and 
other non-forest objects) inside the segments with forest. 
Generally, the tree crown is not spectrally homogenous. 
DN values gradually decrease towards crown edges. 
This is caused by the fact that the top of tree crown is 
more illuminated and crown edges are more shadowed. 
We distinguished tree crowns by applying pixel-based 
approach. We used the method of Isodata that belongs to 
the group of unsupervised classifications. Isodata clas- 
sification was performed independently for each of 13 
photographs. We specified 9 output classes and the min- 
imum cluster size: 100 pixels (fig. 5c). Subsequently, we 
assigned a land cover type to each class. We reclassified 
the result into the values of 1 - forest and 0 - other cat- 
egories. The forest mask was obtained with reciprocal 
combination of outputs from object- and pixel-based ap- 
proaches to the classification of each of 13 photographs. 



First we evaluated the quality of forest masks visually, 
especially comparing them with photographs. Tree 
crown were underestimated on the photographs 10 and 
12. We enlarged the forest mask (surroundings of clas- 
sified tree crown) by about 1 pixel in these photos. For- 
est was overestimated on the photographs 8, 14, 16, 24 
and 26, therefore we reduced the mask by about 1 pixel. 
Overestimation was caused by inclusion of shadowed 
parts of crowns into the mask. The technique of picture 
filtration was used for overestimation elimination with 
kernel size 3 x 3 pixels. No filtering was needed for the 
photographs 4, 6, 18, 20, 22, 28. 

In the next step we created a mosaic from 13 masks. 
In overlap we took the mask with higher quality accord- 
ing to a visual examination. 

Accuracy of the forest mask 

The accuracy of the forest mask was verified on the 
sample of 100 systematically selected points in a grid 
lxl km on the mosaic a subset of segments with for- 
est. Correctness of each point assignment into forest or 
forestless category was visually assessed. The result 
was processed in a form of error matrix. 

Regression model for damage estimation 

The cornerstone for the derivation of tree damage from 
the aerial photograph is the intensity of reflected elec- 
tromagnetic radiation recorded by a digital multispec- 
tral camera. A close relation between defoliation evalu- 
ated in the field survey and the value of spectral reflec- 
tivity recorded with sensors have been proved in several 
works (e.g. Hildebrandt et al. 1991; Gross et al. 2000; 
Rasi and Bucha, 2001; Bucha et al. 2009). Forest health 
evaluation was carried out in two stages. 

First stage 

In the first stage we derived two components using 
orthogonal transformation of original R, G, B and IR 
bands. Mathematic apparatus for the components der- 
ivation was the same as described by Jackson (1983). 
There is a difference in the way of defining the physical 
meaning of the derived components. In Jackson’s ap- 
proach, the first component represents brightness. It was 
derived from the points with low and high reflectivity, 
represented by moist and dry soil. The second compo- 
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nent (greenness) represents the amount of green vegeta- 
tion. In our approach the first component (new synthetic 
channel - further NSCl) emphasizes spectral variability 
of individual tree species. The second component (fur- 
ther NSC2 ) is optimized for damage estimation. 

We chose fully foliaged woods (black poplar and 
willow) with very different reflectivity to compute the 
NSCl and we calculated the differences of DN values: 

b i = (X w -X bp)i 

where: 

i - 1 st to 4 th channel of the aerial photo (IR, R, G, B), 
X - DN value of willow and black poplar. 



Low DN values represent black poplars and Eu- 
roamerican poplars that have low reflectivity in the 
spectrum observed ( B , G, R and IR). The highest DN 
values are represented by willows and white poplars 
which usually indicate high reflectivity. DN values (low 
and high) were calculated as mean values for several 
poplar and willow crowns distinguished manually on 
the aerial photos. 

We standardized the vector b = ( b h b 2 , b 3 , b 4 ) into 
the unit vector by dividing each of his elements by the 
normalization factor B\ 



A u = b/B 



where: 

( 4 

B = 



1? 



\i i y 

A ! ,■ - are sought coefficients for the first component 
derivation: 



NSCl = A U1 x X, + A u2 x X 2 + A X3 *x 3 +A u 4 x X 4 
where X '■ are DN values of pixel in i-th channel. 



While deriving the NSCl we defined DN value rep- 
resenting the maximum damage (dead tree) and calcu- 
lated the difference to any of DN values lying on the 
line between X w and X bp . We chose the value of X w for 
a calculation: 

Si )X dead tree ^w)i ^2,\ * -^1 ,/ 

where D 21 = ^(x^ - X w ) x A iJ 

i 1 



Orthogonality of b and g vectors is ensured with 
described process. Standardization of g vector was car- 
ried out by means of normalisation factor G, and trans- 
formation coefficients were calculated as follows: 



where G = 



( 4 A 

is- 

v -■ i y 



Si/G 



Calculation of the NSC2 was carried out as follows: 
NSC2 = A 2 l xXi + A 2 2 x X 2 + A 2 3 x X 3 + A 2 4 x X 4 

Calculated values of NSCl component represent 
a perpendicular distance from the line of DN values of 
component NSCl defined by the fully foliated trees of 
willow and poplar on its borders. The distance from the 
line is proportional to the extent of the tree damage in 
the given pixel. 



Second stage 

The second stage of forest health evaluation concerned 
obtaining data on tree defoliation during the field sur- 
vey or based on visual interpretation on the monitor 
screen carried out in 3D environment (evaluation of de- 
foliation on stereo-pairs of photographs). For the needs 
of vegetation health classification we evaluated 38 trees 
on stereo-pairs of photographs. 

Predominant and dominant trees were chosen. 
Each tree crown was manually vectorised and mean 
spectral reflectivity was calculated. This spectral char- 
acteristic was interconnected with defoliation data, 
and thus we obtained a set of pair values for further 
analysis. We derived the regression model by means 
of simple linear regression analysis between the data 
obtained both in the first stage ( NSCl component de- 
rived from aerial photographs) and in the second stage 
(defoliation evaluated in the terrain or from stereo 
photographs). Finally, we estimated the damage for 
each picture element (pixel) of NSCl using regression 
equation (Smelko 1990). 



Results 

Forest mask accuracy 

The forest mask accuracy was verified on the sample 
of 100 systematically selected picture elements (tab. 2) 
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and expressed in a form of error matrix (Congalton et 
al. 1999). 



Tab. 2. Error matrix of the forest mask 



Ground truth: 
Classification result 


Forest 


Non- 

forest 


I 


User’s 

Accu- 

racy 


Error 
of com- 
mission 


Forest 


26 


3 


29 


89.7% 


10.3% 


Non-forest 


8 


63 


71 


88.7% 


11.3% 


I 


34 


66 


100 






Producer’s Accuracy 


76.5% 


95.5% 




89.0% 




Error of omission 


23.5% 


4.5% 









Tab. 2 shows that forest is undervalued. In 8 out of 
34 cases, forest was incorrectly incorporated into the 
category of forestless, which represents 23.5% under- 
estimation. 

On the other hand, there were 3 incorrect cases out 
of 29 of forest classifications when forestless was in- 
correctly categorized to forest. This represents 10.3% 
overestimation of category forest 

The overall accuracy of forest mask creation was 
89.0%. We consider this result as acceptable, but not op- 
timal. It points out a need for further refinement of the 
forest mask delineation algorithm from the photographs 
with very high resolution. 

Transformation of aerial photograph bands 
for optimization of damage classification 

Measuring leaf spectral reflectivity in laboratory 
conditions 

In laboratory conditions, with spectral measuring in 
wave spectrum from 450 to 1100 nm there was obtained 
of assimilation apparatus reflectivity of the main tree 
species in the floodplain forest (fig. 3). 

It follows from the results of measurements that 
wavelengths in the red spectrum (630-690 nm) and in 
the infrared spectrum (750 to 770 nm) are suitable for 
the distinction of dead trees from healthy ones. 

Bark reflectivity from dead poplar was higher in the 
blue and green part of the spectrum than the reflectivity 
of willow and poplar leaves’ top side and lower (similar) 
than the reflectivity of the bottom side of willow (pop- 
lar) leaves. This fact can noticeably influence utilization 
of the blue and green channel in damage classification, 



for example in case of photographs taken in windy con- 
ditions. Then, there can be expected higher influence 
of leaf movements on overall reflectivity and spectral 
display of leaves and bark (dead trees) may become un- 
distinguishable in the aerial photograph. 




Fig. 3. Spectral reflectivity of the main tree species leaves and 
tree bark. Top: top of leaf; bottom: bottom of leaf 



Derivation of transformation coefficients from DN 
values of the aerial photograph 

Tab. 3 shows DN values (reflectivity) for willow, poplar, 
and dead trees, based on which we derived coefficients 
for computing NSC 1 and NSC2 components. 



Tab. 3. Input DN values derived from the aerial photograph 





IR 


Red 


Green 


Blue 




DA values for poplar, willow and dead tree 
according to multispectral bands 


Salix alba 


211.2 


116.9 


120.1 


95.6 


Populus nigra 


207.8 


75.8 


81.2 


61.2 


Dead trees 


182.2 


132.4 


123.9 


102.5 


Calculation of transformation coefficient of 1 st and 2 nd 
component 


NSC 1 


0.045 


0.615 


0.585 


0.525 


NSC2 


-0.968 


0.209 


-0.136 


-0.011 



Similar relationships between DN values of the aer- 
ial multispectral bands and reflectivity measured with 
a Licor apparatus were recorded for white willow, black 
poplar and dead tree bark. In the visible part of the spec- 
trum ( R , G and B) reflectivity of Populus nigra < Salix 
alba < dead tree. In the IR channel (755 nm) the relation 
is as follows: dead tree < Populus nigra < Salix alba. 
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Calculation of the components was carried out 
as a linear combination of original aerial photograph 
bands (fig. 5b). We used transformation coefficients for 
the photographs 4 to 24: 

NSC2 a 24 = -0.968 xIR + 0.209 x R + 

-0.136 x G- 0.011 x B [1] 

For the photographs 26 to 28 we used transforma- 
tion as follows: 

NSC2 26 _ 2 z = -0.852 x 1R + 0.064 x R + 

-0.089 x G + 0.512 x B [2] 

where: 

R - red channel, 

G - green channel, 

B - blue channel, 

IR - infrared channel. 

The reason for application of two equations was the 
spectral difference between photographs that is prob- 
ably related to the process of their pre-processing. Dif- 
ferent contribution of the blue channel in equations [1] 
and [2] is mainly visible. 

Derivation of regression model for damage 
estimation 

The results of regression analysis carried out on pair 
values of the first ( NSC2 values) and second stages (de- 
foliation) are regression models (fig. 4), with the help of 
which we estimated damage for each picture element 
(fig. 5d - right bottom) of the first stage: 

Defoliation 4 _ 24 = 292.7 - 1.393 x NSC2 4 _ 14 [3] 

The correctness of the classification can be assessed 
based on the parameters of regression analysis (cor- 
relation coefficient, standard error of regression line). 
At relation [3] the correlation coefficient was r = 0.93 
and standard error s yx = ±13.3%, average deviation 
d t = 10.6%, the range of selection n = 38 trees. Interpre- 
tation of regression line mean error can be illustrated by 
the following example: if derived defoliation is 30%, in 
reality it can deviate in a range ±13.3%, i.e. from 16.7% 
up to 43.3% with 68% reliability. 

When continuous decline rating was rounded to 
the nearest integer for class comparison (5 classes - see 
tab. 4 for class range), the model was able to predict 
decline with 58% accuracy. 




Component NSC2 



Fig. 4. Graph of a regression model [3] between pair value of 
first and second stage of defoliation estimation 

The parameters (coefficients a, b) of linear regres- 
sion in equation [4] for the photographs 26 and 28 was 
deduced in a simplified way, only from trees with no 
defoliation and dead trees (i.e. tree with 0% and 100% 
defoliation). 

Defoliation 26 - 2 i = 386.0 - 2.780 x NSC2 26 _ 2 s [4] 

Due to simplify way of damage estimation a corre- 
lation and mean error of regression line are not derived. 

Classification of forest damage 

Information about damage of forests in the area affect- 
ed by the Gabcikovo waterworks is presented in tab. 4 
and 5. The scale in tab. 4 results from the classification 
used in the ICP Forests Program (UN/ECE ICP Forests, 
2006). Defoliation was assessed for each pixel classified 
as forest in the process of forest mask creation. 

Defoliation calculated for each pixel according to 
equations [3] and [4] was a basis for computation of 
damage in forest compartments. Intensity of damage at 
a level of forest compartment was calculated as arith- 
metic mean of all forest pixels in the compartment. The 
summary results are presented in tab. 5. 

In the area investigated, there were only 2.3% of 
medium damaged stands - mainly in the north-western 
part of the study area, where a decrease of groundwa- 
ter level caused degradation of woody plants. Fleavily 
damaged stands with defoliation above 50% were not 
observed at all. 



Folia Forestalia Polonico, series A, 2013, Vo I. 55 (2), 58-71 



Improved methods of classification of multispectral aerial photographs... 



67 




Fig. 5. A: multispectral aerial photograph, combination of bands: IR/R/G. Cultivated healthy poplar trees are shown in red. Red- 
gray colour with variable texture indicates medium to heavily defoliated tree crowns. Deadwood is shown in colours from gray 
to turquoise. B: NSC2 component derived by orthogonal transformation of original bands of aerial photograph. C: Classification 
of aerial photographs into 9 categories by Isodata method. D: Classification of tree crown defoliation. Green: healthy vegetation. 
Red: damaged vegetation 
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Tab. 4. Summary of the results on classification at a pixel 
level in line with UN-ECE ICP Forests methodology (whole 
observed area) 



Defoliation 

class 


Loss of 
assimilation 
organs (%) 


% of the pixels 
in defoliation 
class 


Description 
of damage 


0 


0-10 


38.0 


No defoliation 


1 


11-25 


37.3 


Slight 

defoliation 


2 


26-60 


23.3 


Moderate 

defoliation 


3 


61-90 


1.3 


Severe 

defoliation 


4 


91-100 


0.1 


Dying and 
dead 


Total 




109 





Tab. 5. Summary results of classification at a forest 
compartment level: frequency of compartments 
in the defoliation classes in the year 2011 



Stand 

defoliation 

(%) 


Description of 
damage 


Number of 
compartments 


Frequency 

distribution 

(%) 


0-10 


Healthy stands 


38 


5.3 


11-20 


Healthy stands 
with first signs of 
damage 


482 


67.1 


21-30 


Slightly damaged 
stands 


181 


25.2 


31-40 


Moderately 
damaged stands 


16 


2.2 


41-50 


1 


0.1 


51-100 


Severely damaged 
stands 


- 


- 




Together 


718 


100.0 




Not evaluated 
stands 


77 


- 



Discussion 

A multitask procedure was applied for large area survey 
of forest condition from multispectral photographs with 
high resolution (30 x 30 cm). Firstly, we focused our re- 
search on accurate distinguishing forest from other land 
cover classes for the reason that all the characteristics 
of forest health were subsequently derived only from 
pixels classified as forest. We suggested the method 



combining object- and pixel-based approaches for the 
following purposes: 

Object-based approach allowed a time efficient 
classification of large non-forest segments with the 
use of combined digital and visual classification. One 
photograph was processed in 100 minutes. The whole 
area from Dobrohosf to Sap (13 photographs), includ- 
ing Hungarian Danube area was classified into forest 
and forestless categories in 22 hours (i.e. 3 working 
days). Further rationalisation is possible by applying 
a fully digital algorithm. However, at average size of 
segment of about 4.8 ha, there can be anticipated a quite 
complex algorithm due to heterogeneity inside the seg- 
ments. Therefore, fast visual classification of the seg- 
ments based on expert knowledge could be considered 
as a practical and operative solution. 

Pixel-based approach was applied as unsupervised 
classification with Isodata methods for tree crown de- 
lineation inside the segments classified as forest. We 
successfully excluded the most of non-forest vegetation 
(crops, meadows, bushes, water, shadows) with overes- 
timation error 10.3% in the category forest. Underesti- 
mation error of forest category equal 23.5% shows a ne- 
cessity for developing an improved algorithm focused in 
particular on tree crown delineation. 

Secondly our research was focused on to the meth- 
od for prediction of forest damage. We proposed the 
procedure based on orthogonal transformation of aerial 
multispectral bands instead of commonly used spectral 
indices (Carter 1993, 1994; Sims and Gamon 2002; Vo- 
gelmann et al. 1993). The purpose of the transforma- 
tion was a reduction of four dimensional data into two 
dimensional space with the aim to optimize informa- 
tion contents. The first component (ASCI) of orthogo- 
nal transformation emphasizes spectral variability of 
individual tree species. The main transformation output 
is NSC2 component optimized for damage estimation. 
Consequently, this is a biophysical indicator, the values 
of which could be recalculated to defoliation through re- 
gression with field assessment of defoliation on sample 
of trees. In the derived regression model concerning 
damage estimation, the value of correlation coefficient 
was r = 0.93 (r-square = 0.86) and the standard error s yx 
= ±13.3%. Validation of the continuous decline predic- 
tion carried out by Pontius et al. (2008) resulted in r- 
square = 0.71 and RMSE = 0.582 for ash forest. Eigirdas 
et al. (2013) achieved the highest correlation coefficients 
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as follows: 0.58 (spruce), 0.60 (pine) and 0.39 (birch) be- 
tween field-estimated and predicted crown defoliations 
from Vexcel UltraCam D photographs with 0.5 x 0.5 m 
resolution. They improved the accuracy by aggregating 
predicted and field estimated tree crown defoliation val- 
ues up to the sample plot level. Correlation coefficients 
between plotwise average values of field-estimated and 
predicted defoliations were around 0.8. 

When continuous decline rating was expressed in 
5 defoliation classes, our model was able to predict de- 
cline for calibration data with 58% accuracy and with 
one-class tolerance with accuracy 95%. In case of 10 
classes with 10% defoliation step, the accuracy would 
be 39%. Pontius et al. (2008) declared 63% accuracy us- 
ing linear regression model with 6 variables (6 chloro- 
phyll and canopy water content sensitive indices) for the 
assessment of ash decline (10 defoliation classes) from 
hyperspectral imagery. 

When compared to the previous monitoring survey 
(Rasi and Bucha 2009) the method proposed is advanta- 
geous because of providing possibility to perform forest 
damage assessment based solely on NSC2 component 
derived from aerial data, excluding the need of addi- 
tional information from the forest management plan. 

The method was also verified at an operational level. 
The study shows that the proposed approach can be used 
to produce detailed maps of flood forest decline embrac- 
ing relatively large areas, such as 3000 ha of forest evalu- 
ated in the Slovak side of the Danube in this study. For- 
est monitoring carried out in 3 -year cycle of repetition 
has proved to be sufficiently accurate and - at the same 
time - economically satisfactory. The total cost of photo- 
graphs of 20 thousand ha area was 4000 €, i.e. 0.2 €/ha. 

Yet, the study showed some problems with process- 
ing aerial imagery. Although all the photographs were 
obtained from one flight line, it was necessary to apply 
two models for /V.S'C'2 derivation. Consequently, there was 
a need to derive two forest decline prediction models, one 
for the photographs 4-24 and the second - for the photo- 
graphs 26-28. The reasons of that could be various, e.g. 
different quality of images pre-processing or different at- 
mospheric conditions during image acquisition in the up- 
per part of the study area. This suggest the need to exam- 
ine an use of original multispectral bands in 90 x 90 cm 
resolution for the classification, instead of pre-processed 
photographs at the level 3, i.e. spatially modified into 30 
x 30 cm resolution according to the panchromatic band. 



From forestry point of view, the continuous mon- 
itoring in the study area disproved the hypothesis on 
large-scale and harmful effects of water construction 
operation on forests situated between Danube original 
bed and supplying, resp. drainage channel. A version of 
hydro-technical measures with damming of arm sys- 
tem and water supplying object with volume of water 
approx. 30 m 3 /s, ensures, in most of the area, suitable 
conditions for healthy development and expected pro- 
duction of wood in floodplain forests. 

Conclusions 

The paper describes the improved procedure of survey- 
ing floodplain forest condition in the Gabcikovo region 
based on the classification of digital multispectral pho- 
tographs. First, we proposed easily practicable method 
combining object- and pixel-based approaches for de- 
lineation of the part of tree crowns for damage evalu- 
ations. Next, the NSC2 component, optimized for the 
damage estimation, was created by orthogonal trans- 
formation of the original bands of multispectral photo- 
graphs. Then, we used the component in conjunction 
with ground control data to predict forest health condi- 
tion over large contiguous areas. The results confirmed 
that NS C2 component was a suitable biophysical indica- 
tor to provide precise information about forest condi- 
tion. Furthermore, forest damage assessment was based 
solely on NSC2 component without using additional in- 
formation from the forest management plan. 
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